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Abstract 


We find that the anomalous dimension (in the MS scheme) of the chromo-magnetic 
operator gghya in the heavy-quark effective theory is 


Tmag — 


Ca^s 


1 + 


-[^CA--^TFn,\-+0(a, 


generalizing the one-loop expression known previously. To derive the two-loop 
result, we use the reparametrization invariance and the virial theorem. After 
performing infrared subtractions, all two-loop integrals are of propagator type and 
are evaluated by a recurrence relation for tensor integrals. 
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1 Introduction 


The heavy-quark effective theory (HQET) is a convenient tool to describe the physics 
of hadrons containing a heavy quark It provides a systematic expansion around 
the limit tuq — oo, in which new symmetries of the strong interactions arise, relating 
the long-distance properties of many observables to a small number of hadronic matrix 
elements. The effective Lagrangian of the HQET is [Q-jQ 

/leff = Ktv-DK + ^ K{iDfK + + 0{l/ml) , (1) 

2mQ Atuq ^ ^ 


where = i[D^, D‘'] is the gluon held-strength tensor, and the velocity-dependent 

held describing a heavy quark inside a hadron moving with velocity v. This held is subject 
to the constraint 'f) hy = h^. The leading term in the ehective Lagrangian, which gives rise 
to the Feynman rules of the HQET, is invariant under a global SU{2nh) spin-havour 
symmetry group, where rih is the number of heavy-quark havours. This symmetry is 
broken by the higher-dimensional operators arising at order I/tuq. The hrst operator 
corresponds to the kinetic energy of the heavy quark inside the hadron, and the second 
operator describes the magnetic interaction of the heavy-quark spin with the gluon field. 
The coefficients Ckm and Cmag result from short-distance effects and, in general, depend 
on the scale at which the operators are renormalized. 

Hadronic matrix elements of the higher-dimensional operators in (|^ play a signihcant 
role in many applications of the HQET. It is therefore important to study the properties 
of these operators, in particular under renormalization. The reparametrization invariance 
of the HQET (an invariance under inhnitesimal changes of the velocity) implies that, 
in regularization schemes with a dimensionless regulator, the kinetic operator is not 
renormalized, i.e. Ckin = 1 to all orders in perturbation theory ||^. (In regularization 
schemes with a dimensionful regulator, the kinetic operator mixes with lower-dimensional 
operators 10-) The renormalization of the chromo-magnetic operator at the one-loop 
order is known for some time [00. It is governed by the anomalous dimension, defined 
as 


bmag 


mag 


^■^mag 

d In/i 


= 2a, 


A. 

da. 


zw 

mag 


( 2 ) 


where Zmag is the renormalization constant that relates the bare operator with the renor¬ 
malized one, Omag = The last relation is valid in dimensional regularization 

with a MS-like subtraction scheme, such as MS 0, and follows from the requirement 
that the anomalous dimension be hnite as e = (4 — (i)/2 —>• 0 0. denotes the 

coefficient of the 1/e pole in Z^ag- At the one-loop order, y^ag = CA(ys/2n 0, 0. In 
this paper, we calculate the anomalous dimension to two-loop order. This completes the 
renormalization of the HQET Lagrangian at order l/mq, and to next-to-leading order 
in as- 


A direct calculation of the renormalization factor Z mag at order would require the 
evaluation of a large number of two-loop diagrams. The equivalent calculation for the 
chromo-magnetic operator of a light quark has been performed in Ref. [ITj. For heavy 
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quarks, we follow a different strategy, which exploits the symmetries of the HQET. 
Instead of the chromo-magnetic operator itself, we study the renormalization of the two 
operators 

or = 9sKT , Of = n^UpOf - n%Of. (3) 

For the special choice T = the hrst operator reduces to the chromo-magnetic 

operator, while the second one vanishes. Because the Feynman rules of the HQFT do 
not involve 7 matrices, the Dirac structure of F is not altered by radiative corrections, 
and it is of advantage to treat it as a general matrix. Then the operators Of and Of 
form a basis of physical (class- 1 ) operators, i.e. operators that do not vanish by the 
equations of motion. The mixing of these operators under renormalization is such that 


n 


or = Zi 07 „. + Z2 07 


2,bare ) 


Of = (Zi + Z 2 


*-^2,bare 


( 4 ) 


In the rest frame of the heavy quark, the operator Of contains only chromo-electric held 
components. The virial theorem of the HQET relates this operator to the kinetic 
operator in (0). Since the kinetic operator is not renormalized, we have the exact relation 
Zi 


Z 2 = 1 , i.e. [11 


Zma.g — Zi — 1 — Z 2 ■ (5) 

Below we calculate Z 2 to order and then use (|^) and (|) to obtain the two-loop 
anomalous dimension of the chromo-magnetic operator. 


2 Two-loop Calculation 

The particular structure of the operator Of greatly facilitates the calculation. To obtain 
the factor Z 2 in (^) at order a^, we calculate the insertion of Of into the amputated 
Green function with two heavy quarks and a (background-held) gluon to two-loop order, 
as well as the one-loop diagrams with counterterm insertions needed to subtract the 
sub divergences of the two-loop diagrams. However, we only need to evaluate those 
diagrams giving a contribution proportional to the tree-level matrix element of Of: 

(Of) = gstaT +vq , (6) 

where v is the heavy-quark velocity, ta a colour matrix, a and a the Lorentz and colour 
indices of the gluon, and q the incoming gluon momentum. It suffices to trace in our 
calculation the terms proportional to v°‘v’^q^, as no other operators yield this Lorentz 
structure. Then, taking into account the Feynman rules of the HQET, it is straightfor¬ 
ward to derive the following selection rules: 

1. If the external gluon is attached to the operator, the diagram is proportional to 
g^la Qj, gua does not contribute to the structure v'^v^q^. 

2. If the external gluon is attached to a heavy-quark line, the diagram only depends 
on n ■ g and does not contribute to the desired structure. On the other hand, the 
external gluon must be connected to a heavy-quark line through other gluons (or 
through a light-fermion or ghost loop), otherwise one cannot get 
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Figure 1: One-loop diagrams contributing to the calculation of the renorma¬ 
lization factor Z 2 . The operator Oi' is represented by a square. 


3. The external gluon must be connected to the operator through other gluons (or 
through a light-fermion or ghost loop). If it is not, one cannot get q^. 

Accounting for these rules, only the two diagrams shown in Fig. need to be calcu¬ 
lated at the one-loop order. They give identical contributions. At the two-loop order, 
all diagrams that yield terms proportional to v°‘v'^q^ can be generated by hrst drawing 
the two-loop diagrams that contribute to the heavy-quark two-point function with an 
insertion of and then attaching an external gluon in all possible ways to one of 
the gluon, ghost or light-quark lines or vertices. This gives rise to the graphs shown in 
Fig. 

In evaluating the two-loop diagrams, we only need to keep terms linear in the gluon 
momentum q. Because the pole parts are polynomial in the external momenta, we can 
hrst take one derivative with respect to q and then set q = 0 and k = fc', so that 
all integrals are of propagator type and depend on the single variable uj = v ■ k = 
V ■ k', where k and k' are the external momenta of the heavy quarks. However, this 
straightforward application of the method of infrared (IR) rearrangement fails for 


some of the diagrams {D 2 , £* 5 , T* 7 , T*i 3 , £* 15 , T*i 8 and D 21 ), for which setting g = 0 after 
differentiation leads to IR divergences. In these cases, we apply a simple variant of the 
so-called R* operation []^ , which compensates these IR poles by a recursive construction 
of counterterms for the IR divergent subgraphs. In our case, after accounting for UV 
and IR counterterms in all subgraphs, the overall counterterm is purely UV and gives 
rise to the desired contribution to the anomalous dimension. 

To illustrate the procedure, consider a typical integral, which arises (as coefficient of 
in the calculation of D 2 '- 


I = / d'^sdH 




{v ■ s + Uj){v ■ t + (x;)(-s 2 )[-(s + - s)^] ' 


(7) 


For q = 0, this integral has an IR divergence when s —0. To construct the IR counter¬ 
term that compensates this singularity, we subtract and add the integral 


/iR= d'^s 


-.2)[-(s-Fg)2] 


dR 




uj{y -tP uj){—py ’ 


( 8 ) 


which is obtained by setting s = 0 whenever possible (i.e. by neglecting s compared with 
t, and V ■ s compared with a;). The two integrals I and 7 ir have the same behaviour 
for s —> 0; their difference is IR hnite and can be evaluated for q = 0. In this case 
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Figure 2: Two-loop diagrams contributing to the calculation of the renormalization 
factor Z^- The notation (2) indicates that a mirror copy of the diagram is included 
implicitly. The shaded circles represent one-loop self-energy insertions. 



Figure 3: Schematic representation of the R* operation. Thick (thin) lines show 
propagators with a large (small) momentum flow. The black dots represent the 
original vertices. 
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JiR becomes a tadpole integral and vanishes, so that (/ — /iR)q=o coincides with the 
original integral with q set to zero. This integral can be calculated using the general 
algorithm described below. On dimensional grounds, it is proportional to (— 2 a;)“^^. 
The contribution -^IR, which is necessary to subtract the IR subdivergence of the original 
integral I, factorizes into an IR counterterm (the s integral) and the original diagram 
with the lines of the IR sensitive subgraph removed. This deletion of lines follows from 
the locality of IR divergences in momentum space. Because of this factorization, IR 
subtraction terms such as Rr are proportional to (— 2 a;)“^'^(—The construction 
just described is schematically shown in Fig. It is necessary that exactly the same 
treatment is followed for the one-loop UV counterterms associated with a given diagram. 
Then the 1 /e poles proportional to ln(— 2 u;//i) and ln(—are cancelled for each 
diagram, and the hnal overall counterterm is local. 

All integrals can be treated in the above manner. The remaining two-loop tensor 
integrals are of the general form 


d'^s dR 


00 


00 






.vs+ 00/ \vt + ooj 


( 9 ) 


The diagrams Dio and D12 can be brought into this form by multiplication with 

1 


1 = 


00 


I (n • s -1- a;) -1- (n ■ f -F a;) - [r> • (s + f) -F a;] I 


( 10 ) 


A standard method is to reduce the tensor integrals to scalar two-loop integrals, which 
can be calculated in an algorithmic way . However, since in our case we need integrals 
with up to hve tensor indices, the reduction to scalar integrals would be rather involved. 
Therefore, we compute the tensor integrals directly. Using the method of integration by 
parts 0 , we obtain the recurrence relation 


(^d — c^i — 0^3 — 2q; 5 “h U/) -|- 0^3 3~^(4 — 5 ) -h oi\ 1^2 

n 

=i: . 

i=i 


({«,}) 




( 11 ) 


which allows us to express any two-loop integral in terms of degenerate integrals, which 
have 02 = 0 , 04 = 0 or 05 = 0 . Here 1 + is an operator raising the index ai by one 
unit etc., and [pj] means that this index is omitted. The degenerate integrals can be 
expressed in terms of products of one-loop tensor integrals, which are given by 

\v ■ s + 00/ ( — s ) 

/ oMl ql-^n 

^ , ( 12 ) 

(-sR [- (s-g)2]^ 
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with 


In{a,P) = 


/3) — 


T{d/2 — a + n) T{2a + (3 — d — n) 

r(a)r(/?) ’ 

V{d/2 -a + n) V{d/2 -p)V{a + (3- d/2) 
V{a)V{p)V{d - a - (3 + n) 


The tensor structures are 

[n/2] 

A'M ■“■(l.ia) = (-l)“-^C(n,j;a) Y.'3 

i=0 Vi=a{ijLi) 

[n/2] 

= X: (-ify D(n,i-,a,ff) Y'9 

i=0 Ui=a(^ii) 


<31^2 _ _ _ gJ32j-lI32j yl'2j + l _ 


V1U2 _ _ _ gy2j-lV2j q’32j + l 


where [n/2] is the largest integer less than or equal to n/2, and 

Ar=A d + 2{n — k — a) ’ 

d + 2{k — 1 — (3) 


C{nJ;a) = H 


D{n,r,a,l3) n t^ci + 2{n-k-a)][d+2{k-a-l3)\ 


( 13 ) 


(14) 


(15) 


The primed sums in (|^ mean a summation over all ^(2j — 1)!! permutations that 
lead to a different assignment of indices. 

Using this technique, we have calculated the pole parts of the two-loop diagrams 
shown in Fig. |^. We adopt the background-field formalism and work in the 

‘t Hooft-Feynman gauge. The final expression for the sum of all diagrams is gauge 
independent. The pole parts proportional to are given in Tab. Q for each diagram. 

Here Ca = N, Cp = {N'^ — l)/(2iV), Tp = 1/2 are the colour factors for an SU{N) 
gauge group, and u/ is the number of light-quark flavours. The renormalization scale 
H is introduced by the replacement of the bare coupling constant with the renormalized 
one through the relation = jl^Zggs, with // = in the MS scheme. 


3 Ultraviolet Counterterms 

The two-loop diagrams in Fig. ^ contain sub divergences, which must be subtracted by 
UV counterterms. We first discuss the one-loop counterterms for the heavy-quark and 
gluon propagators and vertices. The corresponding terms in the Lagrangian are 

k^c.t. = dhhviv ■ dh^-Y (dgAl - 5 hhA Kv- Ahy 

- Sbaa + {dgADB^^A'^^ + {dgAl)A^^^B-^\ , (16) 
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Table 1: Pole parts proportional to ggta^ v^v'^q^ of the two-loop diagrams in units 
of (a^/dTr)^. The diagrams Dq, Dg, D 14 , Diq, ^24 and D 25 vanish or do not yield the 
desired Lorentz structure. 
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where is the conventional gluon held, and the background held. In the Feynman 
gauge, the counterterm coefficients are 


Sh, = Zh — I = 


CXo 


4:716 ’ 




ShhA — ZgZhZ^l^ — 1 — {—Ca + 2 ( 7 / 
Sbaa = ZgZj Za ~ 1 = 5 a ; 


a.. 


47re ’ 


(17) 


where Za is the wave-function renormalization factor for the gluon held, and Z^ = Z~ 


holds for the background held [|T^. Note that the Feynman rule for the counterterm for 


the three-gluon vertex with a background held is the same as that for the usual three- 
gluon vertex, but with a diherent renormalization factor. 

In addition, the two-loop diagrams in Fig. ^contain sub divergences that are removed 
by operator counterterms. To hnd these counterterms, we calculate at the one-loop 
order all insertions of the operator Oi'^ into the amputated Green functions with a non¬ 
negative degree of divergence. They have the held content h^h^B, Ji^hyAA, 

hyhyAB, hruh^BB and hyh^^rjf], where r; is a ghost held. It turns out that there are no 
UV divergences in the Green function containing ghost helds, and that the insertion of 
Oi^ into the heavy-quark two-point function vanishes. The remaining divergences can 
be removed by adding counterterms proportional to the original operators Oi^ and 02 ^ 
, as well as to the operators 


m 


C>r = KT ivDK 

+ K (iv ■ Bt) F - n"g^) K , 

or = KT {[iD^, gi - Q^]) K , 

or = v^VpO^" — v'^VpO^^ , 

or = h.F[ggg"] k, 

or = v'^VpOr-v^VpOr. (18) 


Here is the ‘quantum part’ of the gluon held, dehned by the decomposition A^ = 
B^ -j- Q^. Under a local gauge transformation, the background and the quantum held 
transform as 

B^x) U{x) H^(x) U\x) -1- - [d>^U{x)] U\x ), 

9s 

Q^^{x) U{x) Q^x) U\x ), (19) 


with U{x) G SU{N), so that the gluon held A^ obeys the usual transformation law. The 
operators containing the quantum part of the gluon held are needed as insertions inside 






Table 2: Pole parts proportional to gstaX^ of the one-loop diagrams with counter¬ 

term insertions in units of (as/dvr)^. The diagram Cg does not contribute. 



colour factor 



Cl 

CaCf 

2 

£2 


Ca 

c'i 

5 

3£2 



CaTf rif 

4 

3£2 


C 3 

Cl 

5 

3 e 2 



CaTf nj 

4 

3£2 


Q 

c\ 

5 

3£2 



CaTf % 

4 

3^2 


C 5 

CAiCA - 2 Cf) 

1 

£2 


Ce 

Ca{Ca - 2Cf) 

1 

£2 


C 7 

CaCf 

2 


C 9 

Cl 


T 1 2 
£2 ^ £ 

Sum 

Cl 

11 

3£2 

T + 2 
£2 ^ £ 


CaCf 

4 

£2 



CaTf nj 

4 

3£2 



loop diagrams. Since the quantum held transforms in the same way as the covariant 
derivative, they are all gauge invariant. The (class-2) operator which vanishes by 
the equations of motion, has to be included since the two-loop calculation is performed 
off-shell. The one-loop counterterm coefficients 5n are 

cy cy 

6i = {Ca + 2Cf)-^, 62 = 2S3 = 6^ = -Ca-^, ( 20 ) 

47re 47re 

and S 5 = 6 y = 0. The coefficient 6 ^ is irrelevant for our purposes, since insertions of the 
operator into one-loop diagrams do not contribute to the Lorentz structure v'^v’^q^. 

The one-loop diagrams with counterterm insertions, which are required to cancel the 
sub divergences of the two-loop diagrams in Fig. |^, are shown in Fig. The resulting 
pole parts are listed in Tab. Note that most of the counterterm contributions are pure 
1/e^ poles, since the original diagrams in Fig. do not have 0{e^) terms. When the 
counterterm contributions are added to the results for the two-loop diagrams, all non¬ 
local 1/e divergences proportional to ln(—2a;/p) and ln(—g^/p^) cancel. This is a strong 
check of our calculation. For the coefficient of (a<j/47r)^, indicated by the subscript in 
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Ce 

( 1 ) 




Cg 

( 1 ) 



Figure 4: Counterterm insertion diagrams contributing to the calculation of 
the renormalization factor Z^- 


square brackets, times the tree-level matrix element of we obtain 


17\ 




3e2 9ey 


13\ 


3e2 9ey 


( 21 ) 


To obtain Z 2 , we have to account for the wave-function renormalization of the heavy- 
quark helds. (There is no renormalization of the background held, since ||19|| .) 

This removes the CaCf term. Using then relation (^) leads to 


•^mag 1 T 


A^s 

47re 




The coefficient of the 1/e^ pole in Z ^ag obeys the relation 


7 ( 2 ) _ — Z 

mag, [2] 2 mag. 


hi) 


,[1| (^dg,[il • 


( 22 ) 


(23) 


where /do is the hrst coefficient of the [3 function (see (^ 6 ]) below), and the super- and 
subscripts in round (square) brackets have the same meaning as before. Eq. (p3|) is a 
necessary condition for y^ag to be hnite as e —>• 0 |^. For the anomalous dimension of 
the chromo-magnetic operator in the MS scheme, we now obtain from (0) 


bmag 

This is our main result. 


CAas 

2 tt 


17 


13 


1 + ( Y 3 ~ Tq ) — + 0{ag) 


18 


0.0 


TT 


(24) 


4 Conclusions 

Using a relation between renormalization constants that is a consequence of the re- 
parametrization invariance and the virial theorem, we have found an efficient way to 
calculate the two-loop anomalous dimension of the chromo-magnetic operator in the 


10 































HQET. The calculation involves the evaluation of 25 non-vanishing two-loop diagrams. 
Some of these diagrams have IR divergences that must be regulated by keeping the 
external gluon off-shell. After these divergences are subtracted using a simple variant 
of the R* operation, the external momentum can be set to zero and the resulting two- 
loop propagator-type tensor integrals can be evaluated using a recurrence relation, which 
allows to express them in terms of products of one-loop integrals. 

Our result for the anomalous dimension in is rather simple. In our approach, it 
is obvious from the beginning that the ‘abelian’ colour structures Cp and CpTpUf do 
not appear; because of the selection rules given in section H, the result for the anomalous 
dimension is genuinely non-abelian, proportional to Ca- This would not be obvious if the 
calculation were performed in the standard way by studying the mixing of the chromo- 
magnetic operator with itself. More surprisingly, there is no contribution proportional 
to CaCf in the hnal result, and also all terms proportional to tt^ have disappeared. This 
is in contrast with the two-loop anomalous dimension of current operators, where such 
terms do appear [|T^, 

If, instead of an SU{N) gauge theory, we consider the abelian U{1) gauge theory, the 
right-hand side of (^) vanishes, since then Ca = 0 (and Cp = Tp = 1). In our approach, 
it is easy to see that in the abelian theory the anomalous dimension of the magnetic- 
moment operator actually vanishes to all orders in perturbation theory. According to 
the selection rules, in the absence of gauge-boson self-couplings the external photon can 
be connected to the operator and the heavy-quark lines only through a fermion loop with 
at least four photons attached. Such a diagram is convergent except for subgraphs that 
correspond to charge and held renormalization. Consequently, no operator counterterms 
are needed. 

With our two-loop result for ymag at hand, we can evaluate the Wilson coefficient of 
the chromo-magnetic operator in the effective Lagrangian (|^ to next-to-leading order. 
The result is 

+ ... 

(25) 

where Ci = 2{Ca + Cp) is obtained from one-loop matching [^]. The one- and two-loop 
coefficients of the anomalous dimension and (3 function are 


Craa.gi.'^Q / h) 


' a s{mQ) 


7o/2/3o 


1 + 


Ola 


mq) as 
Cl -I- 


[mq) - asifJ^) 


dvr 


dvr 


'foPi 

2/3o^ 


7o — 2Ca , ~ A —^ CaTp rif , 

Po — -^Ca — -TpHf , 

3d 20 

A = y C-l - y CATpUf - ACpTpnj . ( 26 ) 

As an application, we discuss the mass splitting between the ground-state vector and 
pseudoscalar mesons containing a bottom or a charm quark. To leading order in the 
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(27) 


heavy-quark expansion, one finds |]^ 

^ ^ ml. - ml ^ ^ 0 ( 1 / ) ^ 

Crai,g{mc/fJ.) 

For N = 3 colours and rij = 4 light-quark fiavours (as it is appropriate in the region 
between the bottom- and charm-quark masses) 


R 


/ as{mb) 
\as{mc) 


9/25 



7921 asimc) - asirrib) 
3750 TT 




(28) 


where the non-perturbative parameter Kr accounts for higher-order corrections in the 
heavy-quark expansion. Using as{mb) = 0.22 and asirric) = 0.36, we obtain R ^ 0.76, 
in reasonable agreement with the experimental value i?exp = 0.89 ± 0.01. The next-to- 
leading correction in (^) amounts to minus 10% and deteriorates the comparison with 
the experimental value. However, with A/j 200 MeV required for agreement with i?exp 5 
higher-order corrections in the I/tuq expansion remain moderate in size. 

Acknowledgements: We are indebted to A. Grozin for communicating partial results 
of an independent calculation of ymag based on a different technique. One of us (M.N.) 
would like to thank M. Ciuchini, J. Korner and D. Kreimer for useful discussions. G.A. 
acknowledges a grant from the Generalitat Valenciana, and partial support by DGIGYT 
under grant PB94-0080. He also acknowledges the hospitality of the GERN Theory 
Division, where this research was performed. 


References 


[1] For a review, see: M. Neubert, Phys. Rep. 245, 259 (1994); Int. J. Mod. Phys. All, 
4173 (1996). 

[2] H. Georgi, Phys. Lett. B240, 447 (1990). 

[3] E. Eichten and B. Hill, Phys. Lett. B243, 427 (1990). 

[4] A.F. Falk, B. Grinstein and M.E. Luke, Nucl. Phys. B357, 185 (1991). 

[5] M. Luke and A.V. Manohar, Phys. Lett. B286, 348 (1992). 

[6] L. Maiani, G. Martinelli and G.T. Sachrajda, Nucl. Phys. B368, 281 (1992); G. Mar- 
tinelli and G.T. Sachrajda, Phys. Lett. B354, 423 (1995); Nucl. Phys. B478, 660 
(1996). 

[7] M. Neubert, GERN preprint GERN-TH/96-282 (1996) ||hep-ph/9610471|| , to appear 
in Phys. Lett. B. 

[8] W.A. Bardeen, A.J. Buras, D.W. Duke and T. Muta, Phys. Rev. D18, 3998 (1978). 

[9] E.G. Floratos, D.A. Ross and G.T. Sachrajda, Nucl. Phys. B129, 66 (1977). 

[10] M. Misiak and M. Miinz, Phys. Lett. B344, 308 (1995). 


12 













[11] G. Amoros and M. Neubert, CERN preprint CERN-TH/96-336 (1996) ||hep- 


ph/9612298|| , to appear in Phys. Lett. B. Note that onr definition of the renor¬ 


malization factors differs from the one nsed in this reference. 


[12] M. Nenbert, Phys. Lett. B322, 419 (1994). 

[13] A.A. Vladimirov, Teor. Mat. Phys. 43, 210 (1980). 

[14] K.G. Ghetyrkin and F.V. Tkachov, Phys. Lett. B114, 340 (1982); K.G. Ghetyrkin 
and V.A. Smirnov, Phys. Lett. B144, 419 (1984). 

[15] D.J. Broadhurst and A.G. Grozin, Phys. Lett. B267, 105 (1991). 

[16] K.G. Ghetyrkin and F.V. Tkachov, Nncl. Phys. B192, 159 (1981). 

[17] G. ’t Hooft, in: Fnnctional and Probabilistic Methods in Qnantnm Field Theory, 
Proc. XII Winter School of Theoretical Physics, Karpacz 1975, Acta Univ. Wratisl. 
368, 345 (1976). 

[18] D. Bonlware, Phys. Rev. D23, 389 (1981). 

[19] L.F. Abbott, Nncl. Phys. B185, 189 (1981); Acta Phys. Pol. C13, 33 (1982). 

[20] X. Ji and M.J. Mnsolf, Phys. Lett. B257, 409 (1991). 

[21] G.P. Korchemsky and A.V. Radynshkin, Nncl. Phys. B283, 342 (1987); Phys. Lett. 
B279, 359 (1992); G.P. Korchemsky, Mod. Phys. Lett. A4, 1257 (1989). 

[22] W. Kilian, P. Manakos and T. Mannel, Phys. Rev. D48, 1321 (1993). 


13 







